This project seeks to estimate sport fish harvest and releases of rockfish in Alaska waters by improving on the Howard et al. (2020) methods and expand the time series back to 1977 when the statewide harvest survey (SWHS) was first implemented. This is essentially a Bayesian version of the Howard methods that allows for more appropriate and defensible sharing of information between areas, handles missing data in a more appropriate manor, accurately propagates uncertainty throughout the estimation procedure and replaces the Howard decision tree approach to low sample sizes with a hierarchical model. The methods and results for generating harvest estimates are generally consistent between the Bayesian model and the Howard methods. Harvest estimates are consistent with Howard estimates during contemporary times, but may differ based on more appropriate weighting of SWHS and logbook data, including estimating and correcting bias in the SWHS data.

Table 1. Summary of key improvements in reconstructiing sport fish removals of rockfish using the Bayesian model as compared to the Howard et al. (2020) methods.
Issue Howard Bayes
Time series 1999 - present 1977 - present
Bias in SWHS Not explicitly dealt with. Relies on logbook data and ratios of guided/unguided from SWHS data to estimate unguided releases and harvests. Explicitly estimates bias in SWHS harvest and release estimates based on logbook data.
Species composition of releases Assumes that species composition of releases is equal to that of the harvest, which is not evident in the logbook data. Recognizes different release probabilities by species / species assemblage and estimates it from logbook data and bias corrected SWHS data
Sample size limitations Uses sample size threshholds such that when areas fall below those threshholds values are borrowed from nearby areas. Uses a hierarchichacal modelling approach that shares information between areas in the same region. Thus all data is used, even with small sample sizes. This is a more sound method that avoids assumptions and uses all of the data.
Error propogation Error is propogated when variance estimates are available, but there is uncertainty associated with borrowing values from nearby areas, or the assumption of species compositions being identical in harvest and releases, are not dealt with. By breaking the assumption that species composition is equal between harvests and releases, uncertainty in the release estimates is more reflective of the fishery. Furthermore, the hyerarchichal approach more accurately captures uncertainy within and between areas within a region.
**Figure 1.**- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.

Figure 1.- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.


**Figure 2.**- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).

Figure 2.- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).


Results

**Figure X.**- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Figure X.- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Estimate comparison

Since previous estimates of rockfish harvest have been produced these first 3 graphs will be used to show how the modeled estimates compare to the estimates produced earlier. For total rockfish the estimates are in general agreement although differences are noted. These estimates should be more reliable because they include both SWHS and guide logbook data, handle variance more appropriately, use hierarchical distributions when data is missing, directly consider observation error and are produced using reproducible research.

**Figure 3.**- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 3.**- Total rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


Notes from Adam: When looking at only black rockfish the most significant differences are for the Prince William Sound Inside area. I did not spend a great deal of time tracking this down although it looks like the previous version used bad values for \(P_{(black)ayu}\) for at least unguided anglers. For the moment I would ignore the results for BSIA and SOKO2SAP. I think it is possible to give approximate values for these areas but it will require a little more coding which I have yet to do.

**Figure 4.**- Black rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 4.- Black rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


And black rockfish releases…

**Figure 5.**- Black rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 5.- Black rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 6.**- Yellow rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 6.- Yellow rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 7.**- Yellow rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 7.- Yellow rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 8.**- DSR rockfish (excluding yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 8.- DSR rockfish (excluding yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 9.**- DSR rockfish releases (including yelloweye) 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 9.- DSR rockfish releases (including yelloweye) 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 11.**- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 11.- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 12.**- Slope rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 12.- Slope rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Total Biomass Removal Estimates

**Figure 13.**- Black rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 13.- Black rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.



**Figure 14.**- Yellow rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 14.- Yellow rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

**Figure 15.**- Pelagic rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 15.- Pelagic rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 16.**- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 16.- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 17.**- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 17.- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


Model fit

Logbook residuals

**Figure 18.**- Residuals from logbook harvests.

Figure 18.- Residuals from logbook harvests.


SWHS residuals

**Figure 19.**- Residuals from SWHS harvests.

Figure 19.- Residuals from SWHS harvests.



**Figure 20.**- Residual of SWHS releases.

Figure 20.- Residual of SWHS releases.

Parameter estimates

P(Charter)

These histograms show the posterior distribution of the mean percent of rockfish harvested by the charter fleet.

**Figure 21.**- Mean percent of harvest by charter anglers.

Figure 21.- Mean percent of harvest by charter anglers.


When considered annually we see the percent of rockfish harvested by the charter fleet follows our data fairly well although the model smooths out the changes and we just do not have much information about this ratio. Prior to 2011 the percent charter is confounded with SWHS bias and should be mostly discounted.

**Figure 22.**- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

Figure 22.- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

P(Harvest)

These plots show the fitted logistic line to the proportion of caught rockfish that are harvested. These estimates are used for hindcasting catch estimates based on the harvest data in early years when catch estimates are unavailable.


**Figure 23.**- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 23.- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 24.**- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 24.- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 25.**- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.

Figure 25.- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.


**Figure 26.**- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 26.- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.


**Figure 27.**- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 27.- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

SWHS bias

Figure 23 shows the mean estimate for SWHS bias in harvests and releases. Cook Inlet, North Gulf Coast and North Southeast Inside all look pretty good while most other areas have substantial bias. Prince William Sound Inside has the largest bias. Bias in release estimates is substantial and whereas the SWHS appears to underestimate harvests, it appears to greatly overestimates releases by a factor of 2 or more in most areas as derived from logbook reported releases.

**Figure 28.**- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 28.- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.


Our estimates of SWHS harvest bias track observations fairly well when he have guided harvest estimates. The estimates of release bias in the SWHS data track observed patterns to an extent, but appear to smooth these more volatile disagreements with the logbook data. Adam postulated in his initial start on this that some of this could be the result of the estimates of the proportion guided. This value was not modelled with a trend and thus applies a constant estimate when hindcasting. Data on these relationships could greatly improve this model.

**Figure 29.**- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 29.- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.

P(pelagic)

We model the percentage of pelagic rockfish in the harvest because we have the information for charter anglers (via logbooks) starting in 1998. Other than looking at the model estimates you can use Figure 25 to compare the two data streams for pelagic rockfish harvest. In general they are in agreement with major exceptions in Price William Sound inside, Prince William Sound outside (early in the time series) and South Southeast inside.

**Figure 30.**- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

Figure 30.- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

P(black|pelagic)

Note that in Southeast Alaska we only have composition data starting in 2006. Tania dug up old SE data, but it did not provide any useful data for species apportionment. For the most part, P(black|pelagic) is relatively constant across areas, with the exception of Cook Inlet and NSEI in Southeast AK. It may be worth discussing whether the shifts in those areas is a result of improved or changing species identification rather than actual shift in the species composition of the catch.

**Figure 31.**- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

Figure 31.- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

P(yelloweye|non-pelagic / yelloweye|DSR)

**Figure 32.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 32.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

P(DSR|non-pelagic)

**Figure 33.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

Figure 33.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

P(slope|non-pelagic)

**Figure 34.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 34.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.



P(slope|non-pelagic & non-yellowye) For release estimates

## NULL



Weight Fits

**Figure 35.**- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 35.- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 36.**- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 36.- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 37.**- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 37.- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 38.**- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 38.- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 39.**- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 39.- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


Summary of unconverged parameters:

The proportion of unconverged species composition proportional estimates that are near the bounds of 0 or 1 is 0.3306452

The proportion of unconverged true parameters (defined as Rhat > 1.01) that are also in excess of an Rhat of 1.1 is 0.17 while 0.83 are less than 1.1.

Table Xa. Species proportion derived quantities that are not with 0.05 of the 0 and 1 bounds and are not converged relative to an Rhat of 1.01.

variable

area_n

year

user

mean

median

lower

upper

Rhat

sd

skew

area

p_pelagic

13

1

1

0.504

0.506

0.181

0.817

1.017

0.168

-0.012

NSEI

p_pelagic

13

10

1

0.503

0.504

0.189

0.811

1.017

0.166

-0.006

NSEI

p_pelagic

13

10

2

0.285

0.266

0.077

0.600

1.011

0.137

0.138

NSEI

p_pelagic

13

11

1

0.503

0.504

0.193

0.808

1.013

0.164

-0.006

NSEI

p_pelagic

13

11

2

0.287

0.263

0.073

0.610

1.014

0.142

0.169

NSEI

p_pelagic

13

12

1

0.505

0.511

0.193

0.807

1.018

0.165

-0.036

NSEI

p_pelagic

13

12

2

0.290

0.267

0.076

0.608

1.012

0.143

0.160

NSEI

p_pelagic

13

13

1

0.508

0.508

0.202

0.813

1.011

0.164

0.000

NSEI

p_pelagic

13

13

2

0.286

0.264

0.072

0.620

1.010

0.143

0.154

NSEI

p_pelagic

13

14

1

0.489

0.489

0.189

0.792

1.011

0.163

0.000

NSEI

p_pelagic

13

14

2

0.280

0.263

0.076

0.586

1.011

0.134

0.127

NSEI

p_pelagic

13

15

1

0.454

0.447

0.171

0.787

1.014

0.166

0.042

NSEI

p_pelagic

13

15

2

0.262

0.237

0.071

0.578

1.011

0.133

0.188

NSEI

p_pelagic

13

16

1

0.493

0.488

0.193

0.807

1.011

0.163

0.031

NSEI

p_pelagic

13

16

2

0.284

0.263

0.080

0.602

1.011

0.138

0.152

NSEI

p_pelagic

13

18

2

0.290

0.273

0.078

0.592

1.013

0.138

0.123

NSEI

p_pelagic

13

2

1

0.504

0.503

0.192

0.812

1.021

0.168

0.006

NSEI

p_pelagic

13

2

2

0.283

0.264

0.074

0.595

1.022

0.141

0.135

NSEI

p_pelagic

13

20

2

0.294

0.276

0.084

0.611

1.010

0.141

0.128

NSEI

p_pelagic

13

3

1

0.503

0.508

0.187

0.818

1.019

0.170

-0.029

NSEI

p_pelagic

13

3

2

0.283

0.261

0.071

0.607

1.016

0.144

0.153

NSEI

p_pelagic

13

4

1

0.502

0.506

0.190

0.816

1.012

0.168

-0.024

NSEI

p_pelagic

13

4

2

0.288

0.268

0.072

0.613

1.012

0.143

0.140

NSEI

p_pelagic

13

6

1

0.502

0.501

0.187

0.810

1.013

0.166

0.006

NSEI

p_pelagic

13

6

2

0.283

0.264

0.074

0.603

1.015

0.141

0.135

NSEI

p_pelagic

13

7

1

0.506

0.509

0.190

0.811

1.016

0.168

-0.018

NSEI

p_pelagic

13

7

2

0.285

0.263

0.072

0.615

1.011

0.142

0.154

NSEI

p_pelagic

13

8

1

0.506

0.504

0.196

0.813

1.013

0.164

0.012

NSEI

p_pelagic

13

8

2

0.290

0.268

0.073

0.623

1.011

0.145

0.152

NSEI

p_pelagic

13

9

1

0.504

0.505

0.184

0.808

1.021

0.167

-0.006

NSEI

p_pelagic

13

9

2

0.288

0.264

0.077

0.598

1.012

0.142

0.169

NSEI

p_pelagic

14

11

2

0.367

0.354

0.117

0.698

1.010

0.153

0.085

NSEO

p_yellow

14

12

1

0.894

0.917

0.677

0.983

1.018

0.081

-0.283

NSEO

p_pelagic

14

5

1

0.397

0.384

0.141

0.717

1.011

0.152

0.085

NSEO

p_yellow

15

1

2

0.457

0.449

0.165

0.797

1.012

0.168

0.048

SSEI

p_pelagic

15

15

2

0.173

0.155

0.045

0.410

1.018

0.094

0.191

SSEI

p_pelagic

15

16

2

0.187

0.166

0.046

0.445

1.010

0.102

0.207

SSEI

p_pelagic

15

6

2

0.187

0.167

0.047

0.437

1.016

0.103

0.195

SSEI

p_yellow

2

2

1

0.810

0.830

0.554

0.959

1.010

0.106

-0.189

NG

p_yellow

2

5

2

0.747

0.765

0.441

0.939

1.013

0.131

-0.137

NG

p_pelagic

3

1

1

0.538

0.545

0.215

0.839

1.027

0.167

-0.042

PWSI

p_pelagic

3

1

2

0.200

0.179

0.044

0.468

1.033

0.112

0.187

PWSI

p_pelagic

3

10

1

0.542

0.554

0.209

0.836

1.024

0.166

-0.072

PWSI

p_pelagic

3

10

2

0.204

0.183

0.047

0.489

1.033

0.114

0.183

PWSI

p_pelagic

3

11

1

0.545

0.554

0.214

0.837

1.036

0.167

-0.054

PWSI

p_pelagic

3

11

2

0.206

0.180

0.050

0.490

1.024

0.117

0.223

PWSI

p_pelagic

3

12

1

0.540

0.545

0.220

0.835

1.032

0.164

-0.030

PWSI

p_pelagic

3

12

2

0.205

0.185

0.045

0.488

1.016

0.115

0.174

PWSI

p_pelagic

3

13

1

0.541

0.550

0.208

0.839

1.028

0.169

-0.053

PWSI

p_pelagic

3

13

2

0.204

0.182

0.045

0.485

1.028

0.116

0.190

PWSI

p_pelagic

3

14

1

0.551

0.557

0.226

0.841

1.030

0.165

-0.036

PWSI

p_pelagic

3

14

2

0.215

0.196

0.048

0.506

1.020

0.118

0.161

PWSI

p_pelagic

3

15

1

0.536

0.540

0.206

0.828

1.026

0.163

-0.024

PWSI

p_pelagic

3

15

2

0.192

0.170

0.047

0.455

1.017

0.107

0.206

PWSI

p_pelagic

3

16

1

0.540

0.547

0.223

0.826

1.014

0.163

-0.043

PWSI

p_pelagic

3

16

2

0.201

0.180

0.046

0.476

1.018

0.111

0.189

PWSI

p_pelagic

3

17

1

0.541

0.552

0.221

0.828

1.032

0.165

-0.067

PWSI

p_pelagic

3

17

2

0.200

0.177

0.048

0.479

1.026

0.112

0.205

PWSI

p_pelagic

3

18

1

0.553

0.560

0.227

0.836

1.024

0.164

-0.043

PWSI

p_pelagic

3

18

2

0.209

0.186

0.049

0.496

1.018

0.116

0.198

PWSI

p_pelagic

3

19

1

0.549

0.557

0.226

0.841

1.015

0.162

-0.050

PWSI

p_pelagic

3

19

2

0.209

0.188

0.049

0.475

1.013

0.114

0.184

PWSI

p_pelagic

3

2

1

0.540

0.545

0.206

0.847

1.034

0.169

-0.030

PWSI

p_pelagic

3

2

2

0.204

0.179

0.045

0.487

1.032

0.116

0.215

PWSI

p_pelagic

3

3

1

0.542

0.551

0.206

0.828

1.043

0.165

-0.055

PWSI

p_pelagic

3

3

2

0.200

0.182

0.042

0.479

1.022

0.112

0.161

PWSI

p_pelagic

3

4

1

0.546

0.557

0.215

0.843

1.046

0.167

-0.066

PWSI

p_pelagic

3

4

2

0.204

0.181

0.042

0.485

1.027

0.117

0.197

PWSI

p_pelagic

3

5

1

0.536

0.542

0.209

0.834

1.032

0.167

-0.036

PWSI

p_pelagic

3

5

2

0.202

0.179

0.048

0.465

1.029

0.112

0.205

PWSI

p_pelagic

3

6

1

0.540

0.546

0.211

0.834

1.031

0.168

-0.036

PWSI

p_pelagic

3

6

2

0.204

0.183

0.044

0.488

1.020

0.115

0.183

PWSI

p_pelagic

3

7

1

0.541

0.549

0.213

0.838

1.030

0.167

-0.048

PWSI

p_pelagic

3

7

2

0.205

0.187

0.045

0.481

1.020

0.114

0.158

PWSI

p_pelagic

3

8

1

0.539

0.544

0.203

0.834

1.042

0.167

-0.030

PWSI

p_pelagic

3

8

2

0.206

0.182

0.045

0.486

1.020

0.116

0.207

PWSI

p_pelagic

3

9

1

0.545

0.553

0.208

0.835

1.033

0.166

-0.048

PWSI

p_pelagic

3

9

2

0.201

0.178

0.044

0.475

1.027

0.114

0.202

PWSI

p_yellow

4

30

2

0.652

0.662

0.375

0.881

1.010

0.134

-0.075

PWSO

p_yellow

6

45

2

0.825

0.949

0.001

1.000

1.010

0.277

-0.447

SOKO2SAP

p_black

8

40

1

0.513

0.424

0.339

0.943

1.012

0.206

0.432

afognak

p_pelagic

9

23

1

0.885

0.891

0.772

0.962

1.017

0.050

-0.119

eastside

p_pelagic

9

25

1

0.849

0.852

0.730

0.946

1.012

0.057

-0.052

eastside

Table Xb. True parameters with an Rhat in excess of 1.1

variable

area_n

year/species

mean

median

lower

upper

Rhat

sd

skew

year

area

beta0_pH

10

3

0.808

0.613

-0.772

2.288

3.290

0.897

0.217

northeast

beta1_pH

10

3

2.169

2.588

0.000

5.566

3.328

1.760

-0.238

northeast

beta3_pH

10

3

35.645

35.868

29.190

41.879

1.187

2.783

-0.080

northeast

beta0_pH

11

3

2.049

2.430

0.078

3.131

3.582

0.912

-0.418

CSEO

beta0_pH

11

4

0.693

0.740

-0.338

1.453

1.214

0.467

-0.101

CSEO

beta1_pH

11

3

0.849

0.172

0.000

3.100

4.327

1.048

0.646

CSEO

beta1_pH

11

4

2.047

1.995

1.184

3.199

1.136

0.523

0.099

CSEO

beta3_pH

11

3

38.159

39.620

29.596

42.734

1.447

3.855

-0.379

CSEO

beta3_pH

11

4

40.675

40.626

38.813

42.756

1.110

1.113

0.044

CSEO

beta1_pelagic

11

1.056

1.023

0.555

1.598

1.117

0.268

0.123

CSEO

beta2_yellow

11

-3.389

-1.553

-16.179

-0.394

1.174

4.283

-0.429

CSEO

beta0_pH

12

3

2.616

3.254

-0.156

3.999

4.446

1.310

-0.487

EWYKT

beta1_pH

12

3

1.047

0.020

0.000

4.251

6.046

1.524

0.674

EWYKT

beta3_pH

12

3

38.117

39.340

29.462

42.801

1.613

4.235

-0.289

EWYKT

beta0_pH

13

3

1.718

2.092

-0.054

2.677

4.019

0.824

-0.454

NSEI

beta1_pH

13

3

0.688

0.058

0.000

2.632

5.062

0.929

0.678

NSEI

beta3_pH

13

3

38.248

40.086

29.735

42.696

1.462

3.912

-0.470

NSEI

beta0_pelagic

13

-0.002

0.053

-0.748

0.546

1.102

0.347

-0.159

NSEI

beta1_pelagic

13

1.480

1.358

0.397

3.257

1.210

0.745

0.164

NSEI

beta3_pelagic

13

35.571

35.818

27.903

40.713

1.172

3.315

-0.075

NSEI

beta0_pH

14

3

0.698

0.531

-0.656

2.224

2.674

0.823

0.203

NSEO

beta1_pH

14

3

1.557

1.789

0.000

3.080

2.974

0.972

-0.239

NSEO

beta3_pH

14

3

39.513

40.243

30.958

41.915

1.574

2.526

-0.289

NSEO

beta0_pH

15

3

1.316

1.168

-0.211

2.733

3.299

0.888

0.167

SSEI

beta1_pH

15

3

1.347

1.631

0.000

3.134

3.927

1.074

-0.264

SSEI

beta3_pH

15

3

39.181

40.427

30.013

42.349

1.809

3.197

-0.390

SSEI

beta0_pelagic

15

-0.307

-0.282

-0.820

0.018

1.155

0.196

-0.128

SSEI

beta1_pelagic

15

0.141

0.070

0.011

0.597

1.151

0.161

0.442

SSEI

beta0_yellow

15

-1.179

-1.148

-1.852

-0.777

1.226

0.252

-0.123

SSEI

beta1_yellow

15

0.750

0.697

0.023

2.127

1.206

0.479

0.111

SSEI

beta2_yellow

15

-5.202

-3.008

-19.939

-0.154

1.344

5.467

-0.401

SSEI

beta0_pH

16

3

1.929

1.720

-0.094

4.273

2.077

1.230

0.170

SSEO

beta1_pH

16

3

2.414

2.730

0.000

4.648

2.100

1.413

-0.224

SSEO

beta3_pH

16

3

40.649

41.820

30.691

42.722

1.730

3.061

-0.383

SSEO

beta1_pelagic

16

1.150

1.039

0.440

2.750

1.303

0.543

0.205

SSEO

beta3_pelagic

16

32.940

32.838

27.346

39.862

1.155

3.240

0.031

SSEO

beta1_yellow

16

0.190

0.089

0.011

0.878

1.131

0.277

0.364

SSEO

tau_beta0_pH

2

3

0.591

0.323

0.053

2.704

1.155

1.047

0.256

NG

mu_beta1_pH

2

3

2.259

2.463

0.000

6.933

2.011

2.255

-0.090

NG

tau_beta1_pH

2

1

21.634

11.871

0.347

86.284

1.108

24.164

0.404

NG

tau_beta1_pH

2

3

10.859

1.412

0.042

73.006

1.247

19.394

0.487

NG

tau_beta5_pH

2

3

0.878

0.515

0.058

3.696

1.117

1.624

0.224

NG

beta0_yellow

2

0.513

0.543

-0.127

0.880

1.108

0.243

-0.123

NG

beta1_yellow

2

1.106

1.003

0.502

3.239

1.259

0.562

0.183

NG

beta3_yellow

2

28.974

28.782

23.061

33.720

1.145

2.172

0.088

NG

mu_beta0_pH

3

3

1.733

1.816

-0.111

3.431

2.383

1.002

-0.083

PWSI

tau_beta0_pH

3

3

3.394

0.948

0.115

28.737

1.412

9.431

0.259

PWSI

mu_beta1_pH

3

3

0.932

0.216

0.000

3.366

4.955

1.160

0.617

PWSI

tau_beta1_pH

3

3

5.143

0.222

0.043

41.126

1.620

11.469

0.429

PWSI

beta0_pH

3

1

1.244

1.312

0.313

1.708

1.184

0.340

-0.200

PWSI

beta1_pH

3

1

1.979

1.876

1.282

3.174

1.119

0.484

0.213

PWSI

beta3_pH

3

1

32.574

32.760

27.447

36.224

1.153

2.022

-0.092

PWSI

beta0_pelagic

3

0.174

0.215

-0.798

0.719

1.280

0.348

-0.118

PWSI

beta1_pelagic

3

0.672

0.630

0.035

1.735

1.275

0.391

0.107

PWSI

beta3_yellow

4

28.910

27.867

23.717

35.780

1.104

3.203

0.326

PWSO

beta4_pelagic

5

0.523

0.771

-13.899

13.151

1.173

9.067

-0.027

BSAI

beta0_pH

6

1

-0.700

-0.559

-2.840

0.063

1.211

0.622

-0.227

SOKO2SAP

beta1_pH

6

1

2.326

2.139

1.251

4.498

1.127

0.808

0.232

SOKO2SAP

beta4_pH

6

3

0.914

0.716

0.190

2.469

1.265

2.656

0.075

SOKO2SAP

beta4_pH

6

4

0.986

0.732

0.204

2.694

1.283

4.768

0.053

SOKO2SAP

beta0_pH

8

1

-0.590

-0.531

-1.317

-0.029

1.271

0.349

-0.169

afognak

beta0_pH

8

3

0.414

0.429

-0.198

1.004

1.339

0.330

-0.045

afognak

beta1_pH

8

1

2.801

2.694

1.589

4.581

1.365

0.790

0.135

afognak

beta1_pH

8

3

2.119

2.547

0.000

4.969

3.295

1.667

-0.257

afognak

beta2_pH

8

3

5.493

4.211

-2.918

18.728

1.117

5.478

0.234

afognak

beta3_pH

8

1

37.861

38.368

31.983

41.953

1.792

3.048

-0.166

afognak

beta3_pH

8

3

39.426

41.335

29.141

42.602

2.007

3.866

-0.494

afognak

beta0_pH

9

3

-1.536

-1.398

-3.264

-0.271

1.163

0.818

-0.169

eastside

beta1_pH

9

1

2.350

2.162

1.405

4.490

1.320

0.743

0.253

eastside

beta1_pH

9

3

3.725

3.595

2.240

5.833

1.196

0.971

0.134

eastside

beta2_pH

9

1

0.547

0.414

0.132

1.782

1.278

0.582

0.228

eastside

beta3_pH

9

1

32.552

31.924

27.858

41.533

1.237

3.327

0.189

eastside

tau_beta5_black_kod

0.603

0.279

0.046

2.807

1.105

1.340

0.242

logbc_H

15

28

-0.160

-0.160

-0.406

0.085

1.109

0.125

0.000

SSEI

pG

15

28

0.451

0.451

0.351

0.557

1.199

0.052

0.000

SSEI

lambda_H

5

6.980

0.725

0.024

68.528

1.416

18.563

0.337

BSAI

lambda_H

6

4.731

0.170

0.009

39.422

1.148

11.789

0.387

SOKO2SAP

Table Xc. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.01.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.01

741.00

268.00

1,353.00

399.00

Prop. of total

0.73

0.27

0.77

0.23

Prop. median less than 1

0.19

0.24

0.27

0.43

Prop. median less than 10

0.29

0.57

0.50

0.74

Prop. median less than 50

0.38

0.81

0.65

0.80

Median of the median values

466.31

6.61

10.47

1.92

Mean of the median values

1,358.56

266.58

277.21

152.60

Mean skew of posteriors

0.20

0.29

0.25

0.25

Table Xd. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.1.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.1

16.00

28.00

126.00

82.00

Prop. of total

0.36

0.64

0.61

0.39

Prop. median less than 1

0.25

0.68

0.44

0.45

Prop. median less than 10

0.31

0.93

0.59

0.71

Prop. median less than 50

0.62

1.00

0.88

0.73

Median of the median values

25.38

0.45

3.14

2.15

Mean of the median values

610.27

1.96

110.93

60.80

Mean skew of posteriors

0.17

0.24

0.19

0.19

Table Xe. All parameters and derived quantities that are not converged relative to an Rhat of 1.1.

variable

Freq

beta0_pelagic

2

beta0_pH

13

beta0_yellow

1

beta1_pelagic

5

beta1_pH

14

beta1_yellow

3

beta2_pH

2

beta2_yellow

2

beta3_pelagic

2

beta3_pH

12

beta3_yellow

1

beta4_pelagic

1

beta4_pH

2

H_ay

6

H_ayg

4

H_ayu

6

Hb_ay

5

Hb_ayg

3

Hb_ayu

4

Ho_ay

5

Ho_ayg

3

Ho_ayu

8

Hp_ay

5

Hp_ayg

4

Hp_ayu

5

Hs_ayu

3

Hy_ay

1

Hy_ayg

2

Hy_ayu

3

lambda_H

2

mu_beta0_pH

1

mu_beta1_pH

2

pG

1

pH

42

R_ay

1

R_ayg

15

R_ayu

3

Rb_ay

1

Rb_ayg

13

Rb_ayu

2

Rd_ay

1

Rd_ayu

1

re_pelagic

5

re_pH

127

re_yellow

2

Ro_ay

67

Ro_ayg

62

Ro_ayu

54

Rp_ay

1

Rp_ayg

15

Rp_ayu

2

Rs_ay

18

Rs_ayg

19

Rs_ayu

33

Ry_ay

9

Ry_ayg

12

Ry_ayu

21

sd_pH

1

tau_beta0_pH

2

tau_beta1_pH

2

tau_beta5_pH

1

tau_pH

1

Parameter estimates:

Summary Table of Parameter Estimates
Parameter mean sd Lower_CI Median Upper_CI
mu_bc_H[1] -0.011 0.069 -0.140 -0.012 0.135
mu_bc_H[2] 0.060 0.053 -0.037 0.056 0.173
mu_bc_H[3] -0.311 0.066 -0.436 -0.313 -0.177
mu_bc_H[4] -0.934 0.178 -1.293 -0.932 -0.601
mu_bc_H[5] 1.389 1.212 0.122 1.105 4.419
mu_bc_H[6] -2.352 0.316 -2.951 -2.355 -1.721
mu_bc_H[7] -0.460 0.100 -0.668 -0.457 -0.269
mu_bc_H[8] 0.263 0.356 -0.309 0.230 1.041
mu_bc_H[9] -0.249 0.124 -0.482 -0.252 0.005
mu_bc_H[10] 0.078 0.065 -0.044 0.076 0.212
mu_bc_H[11] -0.125 0.038 -0.200 -0.124 -0.052
mu_bc_H[12] -0.253 0.108 -0.482 -0.249 -0.046
mu_bc_H[13] -0.152 0.070 -0.284 -0.153 -0.014
mu_bc_H[14] -0.305 0.088 -0.487 -0.301 -0.142
mu_bc_H[15] -0.333 0.047 -0.420 -0.334 -0.239
mu_bc_H[16] -0.287 0.359 -0.915 -0.311 0.448
mu_bc_R[1] 1.024 0.124 0.777 1.027 1.263
mu_bc_R[2] 1.012 0.073 0.859 1.014 1.150
mu_bc_R[3] 0.283 0.084 0.115 0.284 0.447
mu_bc_R[4] 0.687 0.195 0.278 0.691 1.055
mu_bc_R[5] 1.667 0.372 0.945 1.658 2.425
mu_bc_R[6] -1.431 0.336 -2.090 -1.428 -0.752
mu_bc_R[7] 0.487 0.173 0.124 0.489 0.821
mu_bc_R[8] 0.645 0.143 0.370 0.644 0.928
mu_bc_R[9] 0.391 0.173 0.026 0.400 0.705
mu_bc_R[10] 0.847 0.110 0.635 0.847 1.067
mu_bc_R[11] 1.096 0.072 0.956 1.095 1.238
mu_bc_R[12] 0.801 0.159 0.494 0.798 1.113
mu_bc_R[13] 1.080 0.091 0.895 1.083 1.254
mu_bc_R[14] 0.866 0.119 0.632 0.866 1.102
mu_bc_R[15] 0.817 0.101 0.614 0.816 1.014
mu_bc_R[16] 1.407 0.103 1.202 1.405 1.609
tau_pH[1] 2.817 0.202 2.440 2.809 3.223
tau_pH[2] 1.516 0.140 1.255 1.508 1.795
tau_pH[3] 1.139 0.128 0.896 1.141 1.382
tau_pH[4] 1.702 0.255 1.262 1.683 2.274
tau_pH[5] 0.468 0.072 0.341 0.463 0.624
beta0_pH[1,1] 0.494 0.221 0.049 0.501 0.915
beta0_pH[2,1] 1.415 0.195 0.980 1.417 1.780
beta0_pH[3,1] 1.244 0.340 0.313 1.312 1.708
beta0_pH[4,1] 1.527 0.243 0.991 1.537 1.982
beta0_pH[5,1] -0.887 0.410 -1.839 -0.834 -0.226
beta0_pH[6,1] -0.700 0.622 -2.840 -0.559 0.063
beta0_pH[7,1] 0.702 0.129 0.455 0.703 0.952
beta0_pH[8,1] -0.590 0.349 -1.317 -0.531 -0.029
beta0_pH[9,1] -0.398 0.370 -1.379 -0.351 0.187
beta0_pH[10,1] 0.141 0.213 -0.277 0.143 0.532
beta0_pH[11,1] 0.069 0.185 -0.316 0.077 0.462
beta0_pH[12,1] 0.507 0.211 0.123 0.495 0.946
beta0_pH[13,1] 0.033 0.152 -0.295 0.043 0.311
beta0_pH[14,1] -0.195 0.213 -0.638 -0.175 0.177
beta0_pH[15,1] 0.011 0.199 -0.385 0.016 0.393
beta0_pH[16,1] -0.052 0.203 -0.509 -0.040 0.298
beta0_pH[1,2] 2.639 0.240 2.166 2.640 3.095
beta0_pH[2,2] 3.027 0.185 2.674 3.024 3.395
beta0_pH[3,2] 3.484 0.200 3.107 3.483 3.877
beta0_pH[4,2] 3.290 0.184 2.946 3.285 3.662
beta0_pH[5,2] 4.512 1.398 2.669 4.247 7.936
beta0_pH[6,2] 3.320 0.266 2.791 3.318 3.853
beta0_pH[7,2] 1.998 0.245 1.507 1.998 2.473
beta0_pH[8,2] 2.891 0.219 2.463 2.890 3.323
beta0_pH[9,2] 3.961 0.302 3.355 3.962 4.546
beta0_pH[10,2] 3.960 0.251 3.476 3.960 4.457
beta0_pH[11,2] 2.299 0.151 2.000 2.299 2.588
beta0_pH[12,2] 1.874 0.219 1.433 1.879 2.299
beta0_pH[13,2] 2.555 0.180 2.211 2.550 2.939
beta0_pH[14,2] 2.009 0.162 1.695 2.008 2.333
beta0_pH[15,2] 2.593 0.186 2.220 2.593 2.969
beta0_pH[16,2] 2.722 0.199 2.346 2.715 3.110
beta0_pH[1,3] 1.123 0.406 0.186 1.153 1.823
beta0_pH[2,3] 2.015 0.232 1.564 2.012 2.470
beta0_pH[3,3] 2.656 0.221 2.222 2.650 3.095
beta0_pH[4,3] 2.993 0.235 2.539 2.990 3.454
beta0_pH[5,3] 2.879 1.660 0.682 2.507 7.268
beta0_pH[6,3] 1.242 0.445 0.470 1.215 2.180
beta0_pH[7,3] 0.483 0.241 0.001 0.482 0.949
beta0_pH[8,3] 0.414 0.330 -0.198 0.429 1.004
beta0_pH[9,3] -1.536 0.818 -3.264 -1.398 -0.271
beta0_pH[10,3] 0.808 0.897 -0.772 0.613 2.288
beta0_pH[11,3] 2.049 0.912 0.078 2.430 3.131
beta0_pH[12,3] 2.616 1.310 -0.156 3.254 3.999
beta0_pH[13,3] 1.718 0.824 -0.054 2.092 2.677
beta0_pH[14,3] 0.698 0.823 -0.656 0.531 2.224
beta0_pH[15,3] 1.316 0.888 -0.211 1.168 2.733
beta0_pH[16,3] 1.929 1.230 -0.094 1.720 4.273
beta0_pH[11,4] 0.693 0.467 -0.338 0.740 1.453
beta0_pH[12,4] 0.714 0.712 -0.697 0.663 2.074
beta0_pH[13,4] 0.445 0.493 -0.583 0.483 1.391
beta0_pH[14,4] -0.009 0.439 -0.916 0.000 0.757
beta0_pH[15,4] 0.324 0.502 -0.816 0.394 1.153
beta0_pH[16,4] 1.052 0.738 -0.314 1.016 2.520
beta0_pH[11,5] 2.944 0.355 2.273 2.945 3.644
beta0_pH[12,5] -0.807 0.551 -1.856 -0.823 0.328
beta0_pH[13,5] 2.659 0.413 1.914 2.632 3.561
beta0_pH[14,5] 3.542 0.484 2.602 3.537 4.512
beta0_pH[15,5] 1.570 0.279 1.062 1.550 2.133
beta0_pH[16,5] 0.531 0.203 0.181 0.517 0.976
beta1_pH[1,1] 2.615 0.448 1.845 2.568 3.580
beta1_pH[2,1] 1.721 0.348 1.096 1.696 2.450
beta1_pH[3,1] 1.979 0.484 1.282 1.876 3.174
beta1_pH[4,1] 2.061 0.454 1.351 2.002 3.171
beta1_pH[5,1] 2.426 0.485 1.640 2.361 3.581
beta1_pH[6,1] 2.326 0.808 1.251 2.139 4.498
beta1_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,1] 2.801 0.790 1.589 2.694 4.581
beta1_pH[9,1] 2.350 0.743 1.405 2.162 4.490
beta1_pH[10,1] 2.378 0.427 1.747 2.311 3.460
beta1_pH[11,1] 3.239 0.279 2.712 3.242 3.800
beta1_pH[12,1] 2.961 0.257 2.438 2.966 3.457
beta1_pH[13,1] 2.648 0.237 2.210 2.643 3.112
beta1_pH[14,1] 3.064 0.285 2.544 3.051 3.671
beta1_pH[15,1] 2.538 0.250 2.045 2.534 3.053
beta1_pH[16,1] 3.096 0.247 2.638 3.090 3.609
beta1_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[1,3] 1.760 1.098 0.029 1.680 4.303
beta1_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,3] 2.119 1.667 0.000 2.547 4.969
beta1_pH[9,3] 3.725 0.971 2.240 3.595 5.833
beta1_pH[10,3] 2.169 1.760 0.000 2.588 5.566
beta1_pH[11,3] 0.849 1.048 0.000 0.172 3.100
beta1_pH[12,3] 1.047 1.524 0.000 0.020 4.251
beta1_pH[13,3] 0.688 0.929 0.000 0.058 2.632
beta1_pH[14,3] 1.557 0.972 0.000 1.789 3.080
beta1_pH[15,3] 1.347 1.074 0.000 1.631 3.134
beta1_pH[16,3] 2.414 1.413 0.000 2.730 4.648
beta1_pH[11,4] 2.047 0.523 1.184 1.995 3.199
beta1_pH[12,4] 3.105 0.719 1.794 3.096 4.576
beta1_pH[13,4] 1.999 0.557 0.988 1.954 3.188
beta1_pH[14,4] 2.275 0.471 1.454 2.256 3.278
beta1_pH[15,4] 2.122 0.561 1.096 2.074 3.363
beta1_pH[16,4] 3.356 0.787 1.830 3.378 4.876
beta1_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,1] 0.903 0.601 0.310 0.725 2.802
beta2_pH[2,1] 1.127 0.814 0.240 0.865 3.337
beta2_pH[3,1] 1.001 0.877 0.153 0.689 3.469
beta2_pH[4,1] 0.723 0.600 0.191 0.525 2.744
beta2_pH[5,1] 1.502 1.349 0.275 1.119 5.110
beta2_pH[6,1] 0.552 0.730 0.119 0.352 2.499
beta2_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,1] 0.384 0.426 0.126 0.277 1.233
beta2_pH[9,1] 0.547 0.582 0.132 0.414 1.782
beta2_pH[10,1] 0.843 0.661 0.234 0.664 2.703
beta2_pH[11,1] 0.774 0.217 0.400 0.752 1.254
beta2_pH[12,1] 1.213 0.524 0.657 1.072 2.768
beta2_pH[13,1] 0.962 0.407 0.492 0.883 2.060
beta2_pH[14,1] 1.020 0.419 0.569 0.923 2.133
beta2_pH[15,1] 0.945 0.374 0.493 0.871 1.852
beta2_pH[16,1] 0.937 0.328 0.523 0.885 1.586
beta2_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,3] 6.665 6.925 -3.252 5.273 22.390
beta2_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,3] 5.493 5.478 -2.918 4.211 18.728
beta2_pH[9,3] 4.519 5.390 0.274 1.545 18.631
beta2_pH[10,3] 4.573 5.714 -3.495 2.576 18.929
beta2_pH[11,3] -7.151 5.414 -20.861 -5.873 -0.565
beta2_pH[12,3] -7.378 5.333 -21.390 -6.125 -0.702
beta2_pH[13,3] -7.177 5.325 -21.196 -5.872 -0.726
beta2_pH[14,3] -7.135 5.325 -21.175 -5.618 -1.029
beta2_pH[15,3] -7.196 5.198 -20.564 -5.870 -0.870
beta2_pH[16,3] -7.318 5.128 -20.336 -5.861 -1.274
beta2_pH[11,4] -2.995 3.028 -11.412 -1.996 -0.284
beta2_pH[12,4] -3.556 3.072 -11.903 -2.626 -0.630
beta2_pH[13,4] -3.264 2.981 -11.161 -2.298 -0.476
beta2_pH[14,4] -3.385 2.818 -11.726 -2.457 -0.791
beta2_pH[15,4] -3.466 3.063 -11.980 -2.396 -0.722
beta2_pH[16,4] -3.581 2.715 -11.516 -2.730 -1.081
beta2_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta3_pH[1,1] 34.756 0.982 32.903 34.708 36.740
beta3_pH[2,1] 32.797 1.336 30.427 32.740 35.767
beta3_pH[3,1] 32.574 2.022 27.447 32.760 36.224
beta3_pH[4,1] 32.717 1.913 29.576 32.546 37.046
beta3_pH[5,1] 27.525 1.032 25.797 27.418 30.309
beta3_pH[6,1] 32.753 3.496 25.232 32.771 40.833
beta3_pH[7,1] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,1] 37.861 3.048 31.983 38.368 41.953
beta3_pH[9,1] 32.552 3.327 27.858 31.924 41.533
beta3_pH[10,1] 32.572 1.257 30.531 32.461 35.849
beta3_pH[11,1] 31.305 0.640 30.157 31.288 32.680
beta3_pH[12,1] 30.476 0.472 29.539 30.479 31.370
beta3_pH[13,1] 32.000 0.668 30.547 31.996 33.276
beta3_pH[14,1] 32.202 0.570 31.189 32.167 33.398
beta3_pH[15,1] 30.907 0.682 29.547 30.892 32.276
beta3_pH[16,1] 30.261 0.492 29.291 30.277 31.167
beta3_pH[1,2] 33.896 5.131 25.422 33.744 42.642
beta3_pH[2,2] 34.127 5.196 25.478 34.171 42.551
beta3_pH[3,2] 30.587 7.192 18.694 30.526 42.419
beta3_pH[4,2] 30.571 7.251 18.638 30.664 42.444
beta3_pH[5,2] 30.299 7.191 18.620 30.085 42.306
beta3_pH[6,2] 30.633 7.213 18.609 30.870 42.290
beta3_pH[7,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[9,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[10,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[11,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[12,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[13,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[14,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[15,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[16,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[1,3] 38.145 3.587 28.852 39.737 42.352
beta3_pH[2,3] 33.947 5.173 25.510 33.916 42.512
beta3_pH[3,3] 30.676 7.191 18.678 30.779 42.400
beta3_pH[4,3] 30.790 7.264 18.801 30.800 42.415
beta3_pH[5,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[6,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[7,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,3] 39.426 3.866 29.141 41.335 42.602
beta3_pH[9,3] 32.088 1.826 28.303 32.773 34.465
beta3_pH[10,3] 35.645 2.783 29.190 35.868 41.879
beta3_pH[11,3] 38.159 3.855 29.596 39.620 42.734
beta3_pH[12,3] 38.117 4.235 29.462 39.340 42.801
beta3_pH[13,3] 38.248 3.912 29.735 40.086 42.696
beta3_pH[14,3] 39.513 2.526 30.958 40.243 41.915
beta3_pH[15,3] 39.181 3.197 30.013 40.427 42.349
beta3_pH[16,3] 40.649 3.061 30.691 41.820 42.722
beta3_pH[11,4] 40.675 1.113 38.813 40.626 42.756
beta3_pH[12,4] 41.774 0.943 39.522 42.135 42.853
beta3_pH[13,4] 41.228 1.276 39.584 41.312 42.794
beta3_pH[14,4] 40.665 0.639 39.282 40.742 41.776
beta3_pH[15,4] 40.892 0.696 39.450 40.877 42.219
beta3_pH[16,4] 42.133 0.461 41.074 42.172 42.883
beta3_pH[11,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[12,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[13,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[14,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[15,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[16,5] 25.000 0.000 25.000 25.000 25.000
beta4_pH[1,1] 0.304 0.052 0.207 0.303 0.409
beta4_pH[2,1] 0.248 0.045 0.162 0.246 0.339
beta4_pH[3,1] 0.224 0.042 0.145 0.223 0.309
beta4_pH[4,1] 0.319 0.042 0.236 0.319 0.403
beta4_pH[5,1] 0.425 0.255 0.084 0.386 0.982
beta4_pH[6,1] 0.769 0.487 0.291 0.715 1.523
beta4_pH[7,1] 0.857 0.323 0.353 0.814 1.610
beta4_pH[8,1] 0.212 0.091 0.060 0.208 0.408
beta4_pH[9,1] 0.255 0.102 0.073 0.248 0.473
beta4_pH[10,1] 0.114 0.045 0.037 0.111 0.211
beta4_pH[11,1] 0.462 0.048 0.372 0.462 0.559
beta4_pH[12,1] 0.553 0.058 0.444 0.552 0.669
beta4_pH[13,1] 0.669 0.061 0.553 0.667 0.789
beta4_pH[14,1] 0.733 0.064 0.609 0.733 0.865
beta4_pH[15,1] 0.744 0.070 0.614 0.740 0.894
beta4_pH[16,1] 1.131 0.078 0.985 1.127 1.293
beta4_pH[1,2] 0.388 0.109 0.192 0.382 0.620
beta4_pH[2,2] 0.541 0.078 0.396 0.538 0.704
beta4_pH[3,2] 0.489 0.071 0.351 0.489 0.634
beta4_pH[4,2] 0.641 0.090 0.473 0.640 0.831
beta4_pH[5,2] 0.735 0.912 0.099 0.522 2.727
beta4_pH[6,2] 0.909 0.854 0.179 0.716 2.813
beta4_pH[7,2] 0.880 0.400 0.355 0.799 1.852
beta4_pH[8,2] 0.879 0.228 0.477 0.863 1.370
beta4_pH[9,2] 1.132 1.552 0.297 0.857 3.604
beta4_pH[10,2] 0.481 0.081 0.323 0.480 0.639
beta4_pH[11,2] 0.939 0.097 0.754 0.936 1.137
beta4_pH[12,2] 0.554 0.117 0.344 0.546 0.805
beta4_pH[13,2] 0.911 0.093 0.743 0.907 1.102
beta4_pH[14,2] 0.890 0.100 0.700 0.887 1.090
beta4_pH[15,2] 0.912 0.087 0.752 0.910 1.089
beta4_pH[16,2] 1.239 0.130 0.998 1.236 1.496
beta4_pH[1,3] 0.398 0.116 0.193 0.392 0.648
beta4_pH[2,3] 0.295 0.087 0.137 0.290 0.478
beta4_pH[3,3] 0.270 0.074 0.138 0.267 0.422
beta4_pH[4,3] 0.294 0.076 0.153 0.291 0.450
beta4_pH[5,3] 0.769 0.967 0.117 0.571 2.672
beta4_pH[6,3] 0.914 2.656 0.190 0.716 2.469
beta4_pH[7,3] 0.890 0.430 0.347 0.823 1.860
beta4_pH[8,3] 0.478 0.221 0.135 0.449 1.000
beta4_pH[9,3] 0.896 0.827 0.177 0.690 2.819
beta4_pH[10,3] 0.213 0.085 0.073 0.206 0.401
beta4_pH[11,3] 0.460 0.083 0.302 0.458 0.631
beta4_pH[12,3] 0.439 0.074 0.300 0.437 0.589
beta4_pH[13,3] 0.927 0.132 0.697 0.920 1.214
beta4_pH[14,3] 0.948 0.144 0.687 0.934 1.266
beta4_pH[15,3] 0.670 0.100 0.492 0.666 0.881
beta4_pH[16,3] 0.622 0.088 0.459 0.619 0.804
beta4_pH[1,4] 0.640 0.279 0.262 0.595 1.313
beta4_pH[2,4] 0.711 0.352 0.270 0.640 1.581
beta4_pH[3,4] 0.754 0.420 0.250 0.670 1.749
beta4_pH[4,4] 0.709 0.341 0.265 0.644 1.609
beta4_pH[5,4] 0.780 0.943 0.120 0.565 2.764
beta4_pH[6,4] 0.986 4.768 0.204 0.732 2.694
beta4_pH[7,4] 0.892 0.404 0.365 0.821 1.859
beta4_pH[8,4] 0.730 0.552 0.168 0.616 2.043
beta4_pH[9,4] 0.787 1.039 0.140 0.571 2.615
beta4_pH[10,4] 0.719 0.627 0.154 0.567 2.033
beta4_pH[11,4] 0.587 0.068 0.458 0.586 0.724
beta4_pH[12,4] 0.470 0.070 0.337 0.468 0.611
beta4_pH[13,4] 0.827 0.108 0.638 0.819 1.063
beta4_pH[14,4] 0.980 0.108 0.780 0.976 1.194
beta4_pH[15,4] 0.660 0.080 0.518 0.654 0.839
beta4_pH[16,4] 0.615 0.065 0.490 0.613 0.748
beta4_pH[1,5] 0.638 0.271 0.257 0.588 1.341
beta4_pH[2,5] 0.708 0.332 0.273 0.650 1.521
beta4_pH[3,5] 0.758 0.420 0.238 0.672 1.819
beta4_pH[4,5] 0.714 0.364 0.264 0.649 1.534
beta4_pH[5,5] 0.773 1.032 0.104 0.566 2.487
beta4_pH[6,5] 0.913 0.755 0.211 0.741 2.682
beta4_pH[7,5] 0.899 0.443 0.347 0.823 1.872
beta4_pH[8,5] 0.739 0.644 0.179 0.614 2.004
beta4_pH[9,5] 0.754 0.868 0.123 0.572 2.533
beta4_pH[10,5] 0.718 0.618 0.135 0.565 2.199
beta4_pH[11,5] 0.629 0.115 0.430 0.621 0.884
beta4_pH[12,5] 0.578 0.186 0.284 0.554 1.002
beta4_pH[13,5] 1.352 0.279 0.904 1.312 1.995
beta4_pH[14,5] 0.886 0.144 0.626 0.879 1.183
beta4_pH[15,5] 1.234 0.298 0.765 1.197 1.932
beta4_pH[16,5] 8.891 5.511 3.097 7.511 23.204
beta0_pelagic[1] 2.221 0.136 1.951 2.220 2.487
beta0_pelagic[2] 1.459 0.135 1.168 1.470 1.698
beta0_pelagic[3] 0.174 0.348 -0.798 0.215 0.719
beta0_pelagic[4] 0.155 0.325 -0.502 0.154 0.804
beta0_pelagic[5] 1.434 0.214 1.020 1.436 1.857
beta0_pelagic[6] 1.352 0.155 1.044 1.350 1.655
beta0_pelagic[7] 1.537 0.150 1.236 1.542 1.830
beta0_pelagic[8] 1.772 0.150 1.483 1.771 2.070
beta0_pelagic[9] 2.716 0.163 2.394 2.714 3.044
beta0_pelagic[10] 2.548 0.134 2.284 2.548 2.815
beta0_pelagic[11] -0.017 0.223 -0.446 -0.008 0.437
beta0_pelagic[12] 1.790 0.140 1.515 1.792 2.050
beta0_pelagic[13] -0.002 0.347 -0.748 0.053 0.546
beta0_pelagic[14] -0.486 0.240 -0.957 -0.478 -0.038
beta0_pelagic[15] -0.307 0.196 -0.820 -0.282 0.018
beta0_pelagic[16] 0.004 0.236 -0.466 0.004 0.470
beta1_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[3] 0.672 0.391 0.035 0.630 1.735
beta1_pelagic[4] 0.976 0.339 0.265 0.994 1.607
beta1_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[11] 1.056 0.268 0.555 1.023 1.598
beta1_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[13] 1.480 0.745 0.397 1.358 3.257
beta1_pelagic[14] 1.875 0.445 1.191 1.815 3.011
beta1_pelagic[15] 0.141 0.161 0.011 0.070 0.597
beta1_pelagic[16] 1.150 0.543 0.440 1.039 2.750
beta2_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[3] 6.972 6.058 0.197 5.481 21.283
beta2_pelagic[4] 6.933 5.886 0.310 5.449 21.349
beta2_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[11] 2.148 1.663 0.341 1.656 6.450
beta2_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[13] 1.231 1.610 0.095 0.442 5.625
beta2_pelagic[14] 1.174 1.151 0.235 0.766 4.661
beta2_pelagic[15] 2.497 1.881 0.108 2.142 6.828
beta2_pelagic[16] 1.671 1.644 0.144 1.052 6.172
beta3_pelagic[1] 29.411 6.670 18.570 29.455 40.407
beta3_pelagic[2] 29.286 6.678 18.491 29.361 40.483
beta3_pelagic[3] 28.158 4.061 20.535 28.414 38.361
beta3_pelagic[4] 25.276 1.974 21.328 25.580 28.891
beta3_pelagic[5] 29.663 6.681 18.531 30.209 40.459
beta3_pelagic[6] 29.647 6.711 18.625 29.857 40.445
beta3_pelagic[7] 29.410 6.598 18.502 29.472 40.455
beta3_pelagic[8] 29.412 6.656 18.503 29.403 40.328
beta3_pelagic[9] 29.285 6.620 18.579 29.217 40.415
beta3_pelagic[10] 29.497 6.600 18.537 29.608 40.531
beta3_pelagic[11] 30.541 1.634 27.658 30.443 33.839
beta3_pelagic[12] 29.394 6.630 18.593 29.367 40.462
beta3_pelagic[13] 35.571 3.315 27.903 35.818 40.713
beta3_pelagic[14] 33.188 1.853 30.588 32.861 38.566
beta3_pelagic[15] 30.033 7.046 18.504 30.436 40.536
beta3_pelagic[16] 32.940 3.240 27.346 32.838 39.862
beta0_yellow[1] -0.492 0.162 -0.811 -0.491 -0.179
beta0_yellow[2] 0.513 0.243 -0.127 0.543 0.880
beta0_yellow[3] -0.140 0.212 -0.585 -0.121 0.228
beta0_yellow[4] 0.979 0.246 0.365 1.016 1.366
beta0_yellow[5] 0.418 0.340 -0.350 0.458 0.967
beta0_yellow[6] 0.817 0.165 0.504 0.817 1.148
beta0_yellow[7] 0.932 0.166 0.610 0.928 1.267
beta0_yellow[8] 0.966 0.168 0.656 0.960 1.299
beta0_yellow[9] 0.556 0.174 0.205 0.560 0.880
beta0_yellow[10] 0.698 0.140 0.419 0.700 0.976
beta0_yellow[11] -0.081 0.190 -0.471 -0.072 0.268
beta0_yellow[12] -1.211 0.234 -1.697 -1.199 -0.796
beta0_yellow[13] -1.330 0.267 -1.973 -1.311 -0.880
beta0_yellow[14] -0.274 0.199 -0.690 -0.267 0.101
beta0_yellow[15] -1.179 0.252 -1.852 -1.148 -0.777
beta0_yellow[16] -0.342 0.209 -0.808 -0.327 0.006
beta1_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[2] 1.106 0.562 0.502 1.003 3.239
beta1_yellow[3] 0.410 0.314 0.013 0.396 1.051
beta1_yellow[4] 1.173 0.533 0.467 1.083 2.590
beta1_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[11] 2.048 0.350 1.460 2.018 2.878
beta1_yellow[12] 0.268 0.479 0.011 0.092 1.566
beta1_yellow[13] 2.316 0.656 1.371 2.191 3.960
beta1_yellow[14] 2.672 0.597 1.679 2.604 4.073
beta1_yellow[15] 0.750 0.479 0.023 0.697 2.127
beta1_yellow[16] 0.190 0.277 0.011 0.089 0.878
beta2_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[2] -6.898 6.000 -21.344 -5.554 -0.095
beta2_yellow[3] -7.567 5.820 -21.264 -6.297 -0.194
beta2_yellow[4] -6.515 6.074 -21.102 -4.796 -0.144
beta2_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[11] -3.389 4.283 -16.179 -1.553 -0.394
beta2_yellow[12] -7.910 5.916 -21.513 -6.750 -0.287
beta2_yellow[13] -2.188 3.550 -13.999 -0.748 -0.211
beta2_yellow[14] -4.586 4.478 -16.428 -2.766 -0.515
beta2_yellow[15] -5.202 5.467 -19.939 -3.008 -0.154
beta2_yellow[16] -7.736 5.984 -22.055 -6.528 -0.116
beta3_yellow[1] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[2] 28.974 2.172 23.061 28.782 33.720
beta3_yellow[3] 31.867 4.494 20.771 32.342 40.692
beta3_yellow[4] 28.910 3.203 23.717 27.867 35.780
beta3_yellow[5] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[6] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[7] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[8] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[9] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[10] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[11] 30.172 0.783 28.788 30.065 31.796
beta3_yellow[12] 35.128 4.276 28.364 34.849 42.635
beta3_yellow[13] 30.314 1.086 28.311 30.224 32.667
beta3_yellow[14] 29.356 0.557 28.339 29.268 30.607
beta3_yellow[15] 31.201 2.726 28.123 30.362 38.115
beta3_yellow[16] 34.464 4.087 28.313 33.865 42.465
beta0_black[1] 0.016 0.149 -0.275 0.017 0.305
beta0_black[2] 1.895 0.140 1.624 1.893 2.166
beta0_black[3] 1.356 0.138 1.086 1.355 1.628
beta0_black[4] 2.491 0.147 2.199 2.492 2.772
beta0_black[5] 4.690 2.224 1.825 4.196 10.309
beta0_black[6] 4.546 1.987 2.137 4.027 9.957
beta0_black[7] 3.332 1.781 1.304 2.877 8.360
beta0_black[8] 0.841 0.220 0.432 0.836 1.293
beta0_black[9] 2.600 0.255 2.098 2.594 3.089
beta0_black[10] 1.348 0.139 1.074 1.346 1.616
beta0_black[11] 2.972 0.147 2.689 2.972 3.261
beta0_black[12] 4.772 0.196 4.394 4.773 5.166
beta0_black[13] 0.214 0.171 0.007 0.175 0.629
beta0_black[14] 2.837 0.161 2.526 2.837 3.142
beta0_black[15] 1.284 0.148 0.989 1.287 1.563
beta0_black[16] 3.766 0.154 3.467 3.770 4.059
beta2_black[1] 0.000 0.000 0.000 0.000 0.000
beta2_black[2] 0.000 0.000 0.000 0.000 0.000
beta2_black[3] 0.000 0.000 0.000 0.000 0.000
beta2_black[4] 0.000 0.000 0.000 0.000 0.000
beta2_black[5] 0.000 0.000 0.000 0.000 0.000
beta2_black[6] 0.000 0.000 0.000 0.000 0.000
beta2_black[7] 0.000 0.000 0.000 0.000 0.000
beta2_black[8] 0.000 0.000 0.000 0.000 0.000
beta2_black[9] 0.000 0.000 0.000 0.000 0.000
beta2_black[10] 0.000 0.000 0.000 0.000 0.000
beta2_black[11] 0.000 0.000 0.000 0.000 0.000
beta2_black[12] 0.000 0.000 0.000 0.000 0.000
beta2_black[13] -2.665 1.674 -6.986 -2.256 -0.608
beta2_black[14] 0.000 0.000 0.000 0.000 0.000
beta2_black[15] 0.000 0.000 0.000 0.000 0.000
beta2_black[16] 0.000 0.000 0.000 0.000 0.000
beta3_black[1] 25.000 0.000 25.000 25.000 25.000
beta3_black[2] 25.000 0.000 25.000 25.000 25.000
beta3_black[3] 25.000 0.000 25.000 25.000 25.000
beta3_black[4] 25.000 0.000 25.000 25.000 25.000
beta3_black[5] 25.000 0.000 25.000 25.000 25.000
beta3_black[6] 25.000 0.000 25.000 25.000 25.000
beta3_black[7] 25.000 0.000 25.000 25.000 25.000
beta3_black[8] 25.000 0.000 25.000 25.000 25.000
beta3_black[9] 25.000 0.000 25.000 25.000 25.000
beta3_black[10] 25.000 0.000 25.000 25.000 25.000
beta3_black[11] 25.000 0.000 25.000 25.000 25.000
beta3_black[12] 25.000 0.000 25.000 25.000 25.000
beta3_black[13] 38.958 0.818 37.071 39.097 40.202
beta3_black[14] 25.000 0.000 25.000 25.000 25.000
beta3_black[15] 25.000 0.000 25.000 25.000 25.000
beta3_black[16] 25.000 0.000 25.000 25.000 25.000
beta4_black[1] -0.299 0.196 -0.684 -0.299 0.082
beta4_black[2] 0.262 0.191 -0.117 0.261 0.634
beta4_black[3] -0.975 0.191 -1.354 -0.973 -0.608
beta4_black[4] 0.484 0.224 0.047 0.482 0.927
beta4_black[5] 0.751 1.599 -1.294 0.452 4.357
beta4_black[6] 0.908 1.605 -0.798 0.561 4.819
beta4_black[7] 0.695 1.635 -1.155 0.418 4.071
beta4_black[8] -0.144 0.311 -0.757 -0.139 0.456
beta4_black[9] 0.990 0.842 -0.169 0.813 3.052
beta4_black[10] 0.117 0.182 -0.248 0.117 0.468
beta4_black[11] -0.604 0.201 -1.002 -0.604 -0.207
beta4_black[12] 0.569 0.319 -0.057 0.563 1.213
beta4_black[13] -1.210 0.200 -1.594 -1.210 -0.813
beta4_black[14] 0.055 0.231 -0.387 0.051 0.516
beta4_black[15] -0.998 0.194 -1.373 -0.999 -0.608
beta4_black[16] -0.493 0.211 -0.892 -0.490 -0.066
beta0_dsr[11] -0.617 0.173 -0.950 -0.621 -0.252
beta0_dsr[12] 1.453 0.184 1.097 1.452 1.816
beta0_dsr[13] 0.021 0.170 -0.313 0.023 0.343
beta0_dsr[14] -0.164 0.176 -0.511 -0.161 0.182
beta0_dsr[15] -0.059 0.175 -0.411 -0.064 0.311
beta0_dsr[16] -0.104 0.175 -0.446 -0.106 0.238
beta1_dsr[11] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[12] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[13] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[14] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[15] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[16] 0.000 0.000 0.000 0.000 0.000
beta4_dsr[11] 1.087 0.231 0.632 1.087 1.529
beta4_dsr[12] 0.736 0.334 0.079 0.733 1.419
beta4_dsr[13] 0.220 0.225 -0.222 0.219 0.657
beta4_dsr[14] -0.047 0.257 -0.533 -0.047 0.448
beta4_dsr[15] 0.786 0.231 0.342 0.788 1.253
beta4_dsr[16] 0.461 0.243 0.000 0.457 0.945
beta0_slope[11] -1.914 0.170 -2.239 -1.918 -1.567
beta0_slope[12] -5.469 0.317 -6.119 -5.454 -4.891
beta0_slope[13] -1.203 0.169 -1.547 -1.200 -0.872
beta0_slope[14] -2.565 0.184 -2.925 -2.565 -2.208
beta0_slope[15] -1.204 0.172 -1.527 -1.205 -0.874
beta0_slope[16] -2.304 0.175 -2.628 -2.311 -1.959
beta1_slope[11] 0.000 0.000 0.000 0.000 0.000
beta1_slope[12] 0.000 0.000 0.000 0.000 0.000
beta1_slope[13] 0.000 0.000 0.000 0.000 0.000
beta1_slope[14] 0.000 0.000 0.000 0.000 0.000
beta1_slope[15] 0.000 0.000 0.000 0.000 0.000
beta1_slope[16] 0.000 0.000 0.000 0.000 0.000
beta4_slope[11] -0.775 0.231 -1.214 -0.783 -0.311
beta4_slope[12] -4.253 0.883 -6.091 -4.215 -2.608
beta4_slope[13] -0.084 0.228 -0.533 -0.084 0.361
beta4_slope[14] -0.359 0.276 -0.903 -0.358 0.173
beta4_slope[15] -0.964 0.226 -1.406 -0.965 -0.525
beta4_slope[16] -0.649 0.251 -1.143 -0.650 -0.167
beta5_slope[11] 2.614 0.266 2.103 2.609 3.138
beta5_slope[12] 1.831 0.471 0.954 1.816 2.748
beta5_slope[13] 1.387 0.268 0.858 1.385 1.913
beta5_slope[14] 3.370 0.408 2.590 3.365 4.203
beta5_slope[15] 1.663 0.261 1.147 1.664 2.177
beta5_slope[16] 2.741 0.346 2.088 2.736 3.430
sigma_H[1] 0.207 0.049 0.116 0.204 0.311
sigma_H[2] 0.175 0.028 0.129 0.172 0.236
sigma_H[3] 0.178 0.040 0.107 0.175 0.265
sigma_H[4] 0.456 0.079 0.330 0.446 0.637
sigma_H[5] 0.962 0.204 0.591 0.949 1.389
sigma_H[6] 0.537 0.192 0.170 0.526 0.934
sigma_H[7] 0.312 0.064 0.218 0.303 0.462
sigma_H[8] 0.426 0.092 0.299 0.413 0.630
sigma_H[9] 0.517 0.115 0.336 0.503 0.777
sigma_H[10] 0.213 0.040 0.145 0.209 0.302
sigma_H[11] 0.266 0.044 0.194 0.262 0.364
sigma_H[12] 0.380 0.144 0.200 0.336 0.708
sigma_H[13] 0.220 0.036 0.158 0.216 0.297
sigma_H[14] 0.503 0.089 0.346 0.500 0.696
sigma_H[15] 0.240 0.038 0.176 0.237 0.325
sigma_H[16] 0.220 0.041 0.152 0.215 0.315
lambda_H[1] 1.996 3.057 0.076 0.965 10.109
lambda_H[2] 5.349 4.935 0.489 3.872 19.151
lambda_H[3] 4.038 5.435 0.246 2.265 18.365
lambda_H[4] 0.006 0.005 0.001 0.005 0.018
lambda_H[5] 6.980 18.563 0.024 0.725 68.528
lambda_H[6] 4.731 11.789 0.009 0.170 39.422
lambda_H[7] 0.012 0.008 0.002 0.010 0.032
lambda_H[8] 8.282 11.354 0.053 4.238 38.821
lambda_H[9] 0.016 0.010 0.003 0.013 0.042
lambda_H[10] 0.222 0.271 0.028 0.156 0.755
lambda_H[11] 0.226 0.325 0.011 0.103 1.079
lambda_H[12] 3.581 4.602 0.160 2.115 15.833
lambda_H[13] 3.199 2.796 0.275 2.389 10.860
lambda_H[14] 3.408 4.196 0.228 2.110 14.622
lambda_H[15] 0.020 0.021 0.003 0.014 0.067
lambda_H[16] 0.999 1.416 0.056 0.545 4.872
mu_lambda_H[1] 3.982 1.859 1.036 3.773 8.111
mu_lambda_H[2] 3.718 1.985 0.578 3.569 8.111
mu_lambda_H[3] 3.396 1.840 0.760 3.128 7.542
sigma_lambda_H[1] 7.903 4.307 1.717 7.097 18.095
sigma_lambda_H[2] 8.346 4.828 0.910 7.772 18.383
sigma_lambda_H[3] 6.131 4.030 0.946 5.230 16.260
beta_H[1,1] 6.454 1.276 3.439 6.668 8.413
beta_H[2,1] 9.759 0.609 8.444 9.798 10.892
beta_H[3,1] 7.808 0.837 5.901 7.914 9.184
beta_H[4,1] 10.788 7.468 -4.085 10.918 25.727
beta_H[5,1] -0.122 2.644 -5.505 0.063 4.159
beta_H[6,1] 2.811 4.367 -7.544 4.173 8.282
beta_H[7,1] 0.854 5.702 -11.291 1.150 10.977
beta_H[8,1] 1.443 4.100 -2.523 1.302 3.711
beta_H[9,1] 13.614 5.699 2.349 13.549 25.193
beta_H[10,1] 6.750 1.843 2.902 6.819 10.332
beta_H[11,1] 4.831 3.512 -2.854 5.466 9.779
beta_H[12,1] 2.625 1.091 0.655 2.534 5.012
beta_H[13,1] 9.064 0.868 7.227 9.115 10.512
beta_H[14,1] 2.122 1.036 0.138 2.120 4.119
beta_H[15,1] -6.779 3.506 -13.283 -6.887 0.788
beta_H[16,1] 3.118 2.320 -0.803 2.851 8.575
beta_H[1,2] 7.765 0.268 7.221 7.767 8.272
beta_H[2,2] 9.893 0.160 9.573 9.894 10.200
beta_H[3,2] 8.833 0.199 8.434 8.828 9.221
beta_H[4,2] 2.928 1.421 0.219 2.902 5.812
beta_H[5,2] 1.832 1.047 -0.184 1.828 3.919
beta_H[6,2] 5.911 1.094 3.526 6.041 7.771
beta_H[7,2] 2.492 1.067 0.614 2.418 4.776
beta_H[8,2] 2.970 1.169 1.085 3.145 4.271
beta_H[9,2] 3.153 1.088 0.999 3.123 5.277
beta_H[10,2] 8.011 0.367 7.286 8.023 8.720
beta_H[11,2] 9.827 0.637 8.840 9.715 11.215
beta_H[12,2] 3.948 0.366 3.240 3.940 4.696
beta_H[13,2] 9.135 0.237 8.695 9.127 9.611
beta_H[14,2] 4.023 0.347 3.374 4.016 4.726
beta_H[15,2] 11.502 0.638 10.160 11.522 12.705
beta_H[16,2] 4.547 0.767 3.029 4.539 6.053
beta_H[1,3] 8.455 0.267 7.965 8.443 9.000
beta_H[2,3] 9.885 0.145 9.591 9.886 10.177
beta_H[3,3] 9.523 0.163 9.213 9.522 9.852
beta_H[4,3] -2.089 0.852 -3.839 -2.064 -0.471
beta_H[5,3] 3.889 0.608 2.622 3.894 5.056
beta_H[6,3] 8.387 1.174 6.674 8.137 10.787
beta_H[7,3] -2.616 0.635 -3.909 -2.608 -1.425
beta_H[8,3] 5.308 0.567 4.677 5.211 6.558
beta_H[9,3] -2.667 0.727 -4.106 -2.657 -1.269
beta_H[10,3] 8.553 0.279 8.012 8.551 9.118
beta_H[11,3] 8.499 0.290 7.877 8.530 9.020
beta_H[12,3] 5.268 0.324 4.495 5.307 5.791
beta_H[13,3] 8.873 0.171 8.529 8.875 9.201
beta_H[14,3] 5.774 0.273 5.153 5.797 6.241
beta_H[15,3] 10.295 0.309 9.708 10.287 10.912
beta_H[16,3] 6.506 0.521 5.410 6.544 7.400
beta_H[1,4] 8.075 0.184 7.699 8.084 8.419
beta_H[2,4] 9.876 0.130 9.609 9.883 10.105
beta_H[3,4] 10.020 0.157 9.695 10.031 10.303
beta_H[4,4] 11.813 0.466 10.875 11.813 12.732
beta_H[5,4] 5.548 0.758 4.315 5.459 7.244
beta_H[6,4] 7.067 0.941 5.085 7.279 8.428
beta_H[7,4] 8.566 0.349 7.877 8.575 9.250
beta_H[8,4] 6.696 0.265 6.117 6.722 7.087
beta_H[9,4] 7.414 0.454 6.520 7.422 8.285
beta_H[10,4] 7.500 0.221 7.075 7.492 7.963
beta_H[11,4] 9.509 0.193 9.128 9.509 9.874
beta_H[12,4] 7.258 0.200 6.885 7.253 7.684
beta_H[13,4] 9.059 0.138 8.767 9.061 9.324
beta_H[14,4] 7.824 0.217 7.431 7.819 8.275
beta_H[15,4] 9.449 0.223 8.988 9.449 9.875
beta_H[16,4] 9.316 0.217 8.937 9.303 9.772
beta_H[1,5] 8.876 0.156 8.548 8.880 9.174
beta_H[2,5] 10.680 0.101 10.490 10.675 10.885
beta_H[3,5] 10.833 0.154 10.545 10.828 11.146
beta_H[4,5] 8.339 0.493 7.401 8.344 9.351
beta_H[5,5] 5.175 0.619 3.676 5.257 6.174
beta_H[6,5] 8.998 0.633 8.035 8.890 10.381
beta_H[7,5] 6.577 0.355 5.901 6.571 7.267
beta_H[8,5] 8.247 0.212 7.901 8.227 8.727
beta_H[9,5] 8.200 0.466 7.275 8.193 9.110
beta_H[10,5] 9.996 0.224 9.563 9.997 10.435
beta_H[11,5] 11.596 0.216 11.169 11.595 12.013
beta_H[12,5] 8.615 0.187 8.245 8.612 8.996
beta_H[13,5] 10.119 0.129 9.877 10.118 10.384
beta_H[14,5] 9.247 0.215 8.840 9.236 9.706
beta_H[15,5] 11.336 0.239 10.878 11.334 11.807
beta_H[16,5] 9.897 0.167 9.534 9.906 10.207
beta_H[1,6] 10.146 0.241 9.760 10.109 10.738
beta_H[2,6] 11.467 0.117 11.236 11.468 11.705
beta_H[3,6] 10.561 0.158 10.231 10.565 10.850
beta_H[4,6] 12.925 0.870 11.175 12.939 14.668
beta_H[5,6] 5.689 0.622 4.613 5.654 6.918
beta_H[6,6] 8.688 0.587 7.438 8.725 9.760
beta_H[7,6] 10.565 0.611 9.404 10.566 11.763
beta_H[8,6] 9.512 0.276 8.937 9.530 9.965
beta_H[9,6] 8.253 0.757 6.809 8.252 9.751
beta_H[10,6] 9.139 0.337 8.396 9.159 9.739
beta_H[11,6] 10.695 0.356 9.925 10.723 11.317
beta_H[12,6] 9.376 0.243 8.926 9.366 9.880
beta_H[13,6] 11.066 0.157 10.774 11.065 11.379
beta_H[14,6] 9.932 0.283 9.368 9.936 10.494
beta_H[15,6] 10.515 0.429 9.676 10.519 11.367
beta_H[16,6] 10.622 0.225 10.132 10.633 11.027
beta_H[1,7] 10.162 1.275 6.957 10.483 11.929
beta_H[2,7] 12.118 0.514 11.059 12.144 13.081
beta_H[3,7] 10.092 0.683 8.596 10.143 11.264
beta_H[4,7] 1.344 4.487 -7.655 1.351 10.319
beta_H[5,7] 5.851 2.013 1.830 5.953 9.666
beta_H[6,7] 7.014 3.234 -0.866 8.029 11.246
beta_H[7,7] 7.684 3.105 1.573 7.707 13.581
beta_H[8,7] 10.881 1.044 9.104 10.869 12.756
beta_H[9,7] 5.105 3.889 -2.784 5.202 12.624
beta_H[10,7] 9.927 1.565 7.210 9.831 13.309
beta_H[11,7] 11.371 1.797 8.234 11.182 15.307
beta_H[12,7] 9.757 0.989 7.570 9.843 11.412
beta_H[13,7] 11.678 0.694 10.127 11.731 12.846
beta_H[14,7] 10.628 0.899 8.808 10.673 12.281
beta_H[15,7] 14.045 2.280 9.397 14.009 18.490
beta_H[16,7] 12.521 1.229 10.605 12.330 15.448
beta0_H[1] 8.408 18.001 -28.942 8.724 45.512
beta0_H[2] 10.637 7.723 -5.352 10.625 27.031
beta0_H[3] 9.572 11.248 -13.067 9.760 32.558
beta0_H[4] 11.446 186.541 -359.232 9.557 404.909
beta0_H[5] 4.158 30.913 -59.564 4.490 59.565
beta0_H[6] 6.753 49.712 -104.753 7.625 115.949
beta0_H[7] 7.589 140.905 -280.684 8.234 300.181
beta0_H[8] 6.227 29.445 -24.448 6.622 31.366
beta0_H[9] 4.956 118.326 -227.950 4.886 231.963
beta0_H[10] 8.842 35.933 -66.637 8.862 82.042
beta0_H[11] 8.137 54.673 -107.650 9.724 118.563
beta0_H[12] 6.750 13.186 -19.358 6.851 33.318
beta0_H[13] 10.000 10.569 -11.092 9.834 31.098
beta0_H[14] 7.425 11.844 -16.334 7.292 31.159
beta0_H[15] 6.090 115.837 -234.573 9.233 237.378
beta0_H[16] 8.507 23.065 -39.287 8.617 57.662